Universal measurement of quantum correlations of radiation 
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A measurement technique is proposed which, in principle, allows one to observe the general space- 
time correlation properties of a quantized radiation field. Our method, called balanced homodyne 
correlation measurement, unifies the advantages of balanced homodyne detection with those of 
homodyne correlation measurements. 
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The pioneering photon correlation experiments per- 
formed by Hanbury Brown and Twiss ml in the 1950th 
have stimulated a series of investigations of the correla- 
tion properties of radiation fields. The quantum field the- 
oretical description of the coherence properties of radia- 
tion was introduced by Glauber |^ . In this approach the 
considered correlation functions contain equal powers of 
negative and positive frequency parts of field operators, 
which is closely related to the possibilities of observation 
by photon correlation measurements. 

A complete description of the quantum statistical 
properties of radiation also requires the consideration 
of more general space-time dependent correlation func- 
tions 3, 4, 5, 6], which are composed of unequal powers 
of photon annihilation and creation operators. Phase- 
sensitive correlations of such a type are not directly ac- 
cessible by photoelectric detection. In principle, they 
could be observed by transmitting the radiation, prior 
to the measurement by photodetectors, through an ap- 
propriately chosen nonlinear medium 6] . The realization 
of such measurements is difficult, in particular, when mi- 
croscopic fields to be measured are too weak for causing 
the needed nonlinear effects. 

Usually phase-sensitive radiation properties are mea- 
sured by homodyne detection where the fields to 
be measured are superimposed with a coherent reference 
field, the so-called local oscillator. The methods of bal- 
anced homodyne tomography 's'l and of unbalanced ho- 
modyning allow one to reconstruct the quantum state 
of an effective single-mode radiation field. The recon- 
struction of moments has also been considered lOj. More 
complex methods of multiport homodyningallow one to 
observe space-time dependent correlations [ll|. However, 
involved reconstruction methods are needed and one ob- 
tains only insight in smoothed quantum states, due to 
noise effects caused by imperfect detection, for details 
see ^2 and references therein. 

Some special radiation properties composed of unequal 
orders of annihilation and creation operators can be mea- 
sured by homodyne correlation techniques 0. In this 
case imperfect detection does not contaminate the de- 
tected correlation functions. The method has been fur- 
ther developed for the measurement of arbitrary mo- 
ments of a single- mode radiation field with the aim 



to characterize the nonclassical properties of radiation. 
However, the use of a weak local oscillator makes it dif- 
ficult to determine moments of high orders. 

Why is the accurate determination of space-time de- 
pendent correlation functions of radiation fields, includ- 
ing those composed of unequal powers of annihilation 
and creation operators, of interest? It would lead to new 
possibilities to study the general (high-order) quantum 
coherence properties of radiation sources, including the 
dynamical properties and the spatial irradiation charac- 
teristics. New types of time-dependent nonclassical cor- 
relation properties could be investigated, which general- 
ize known effects like photon antibunching |l5l | . Last but 
not least, the characterization of entanglement of contin- 
uous quantum states can been based on such correlation 
functions 0|. Hence their measurement is of great in- 
terest for any kind of application of nonclassical and, in 
particular, of entangled radiation fields. 

In this letter we propose a method for observing 
the most general normally- and time-ordered correlation 
functions of radiation fields. It combines advantages of 
balanced homodyning with those of homodyne correla- 
tion measurements in a method to be called balanced 
homodyne correlation (BHC) measurement. A chosen 
correlation function is determined by a fixed number of 
photodetectors, so that imperfect detection does not lead 
to smoothing effects. Since a strong local oscillator can 
be used in the BHC method, the signal-to-noise ratio al- 
lows to determine high-order correlation functions. 

Let us consider the general correlation function 
tj(n,m)^-j;i, . . . , Xn, Hi, ... , i/m) of an electromagnetic field. 
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where Xk = (rfc,tfc), yi = (si,r;) refer to both space and 
time points. The operator £^^^ i^^'^'^) denotes the nega- 
tive (positive) frequency part of the electric field operator 
containing the photon creation (annihilation) operators. 
The notation o S , as used in 0| , means that field opera- 
tors are to be written in normal order {£^^^ to the left of 
£ '■"'"•'), and time order (time arguments increasing to the 
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FIG. 1: Basic measurement device MDd. 
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FIG. 2: The recursive structure of MD^ 
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right in products of £^~^ and to the left in products of 
For simphcity we restrict ourselves to the case of 
one polarization, the extension to different polarizations 
is straightforward. 

Let us start with the simplest case of one and the same 
space-time point in all the field operators in the expres- 
sion 1^: xi — . . . = Xn = yi = ■ ■ ■ = Urn = [Y,t). In such 
a case the correlation function reads as 
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Here the field operators are already normally ordered 
and time-ordering becomes meaningless. The correlation 
function ^ can be measured with the device shown in 
Fig. n This device is parameterized by an integer, the 
number of levels or depth of the device. The measure- 
ment device (MD) of the depth d we denote by MD^. It 
has 2'^ photodetectors and can be composed recursively 
of two devices MD^_]^ and MD^'_]^ of lower depth d — 1, 
cf. Fig. 121 The elementary building block of our setup 
is the lowest order device MDi. Roughly speaking, MD^ 
can be constructed by replacing both photodetectors of 
MDi by MDd-i- All the beamsplitters of the device MD^ 
are assumed to be symmetric 50%-50% and all the pho- 
todetectors to have the same quantum efficiency ry. This 
can be realized by balancing them with polarizers. 

In the BHC approach the device MD^ allows us to 
measure the correlations ^ with n + ra ^ 2'^"^, so that 
the minimal depth d necessary to measure the moments 
(0) with given n and m is d = [log2(n -I- to)] -I- 1 ([a;] is 
the smallest integer number greater or equal to x). We 
assume the local oscillator to be in a coherent state 

(3) 



tional to the operators {£ 



, c(±) ^c(±)v 



r^E'+' + vac, E) 



+ vac, 
(4) 

j = 1,...,2'^ ^, the phase depends on the path of 

the signal to the j-th photodetector. The operators fj.^"* 
read as 
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with $_|_ — $_ = tt/2. The symbol vac in Q means 
vacuum terms which play no role in the considered ho- 
modyne correlation measurements. 

Let us denote by ^^'^ normally-ordered (sym- 

bolized by ::) correlation function of the photodetectors 
PD,,,...,PD,,, 
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For any k = 1, . . . , 2'^~^ and for any I, ^ I ^ k, let 
us select I photodetectors from MD'^_]^ and k — I from 
MD'^_i (cf. Fig. 121) and measure their correlation func- 
tion. Due to the relations Q such a correlation function 
depends only on the numbers k and I but not on the in- 

(k) 

dividual photodetectors chosen. We denote it by , 
using Eq. Q it reads as 
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The operators J\f+ and A/"- are defined analogously to the 
photon number operator of a single-mode field, 

Af± = 4-^4+^ = i ± EX^ + E^) , (8) 



One can easily see that the field operators £j (£j ) and corresponds to the quadrature operator. 



detected by j-th photodetector PD^- (PD") of the left 
subdevice MD'^_^ (right subdevice MD^'_j^) are propor- 
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FIG. 3: Measurement scheme for space-time correlations. 

with E^^^ — 5(±)e='=*'^* being slowly varying fields and 
<y5 = "^LO + 7r/2. 

Let us chose one correlation function for each I = 
0, 1, . . . , fc and combine them as 

p^w=^(_i)^-^Qr«. (10) 

It is important to note that all the terms in this sum 
arc proportional to fc-th power of the quantum efficiency 
of the photodetectors. Using the binomial formula we 
obtain 

(11) 

which represents the set of BHC data to be analyzed. The 
dependence of the quantity i^^"') on the phase ip is given 
explicitly: i^*^*^-* — F^'^^ip). According the definition of 
the operator X^p, the normally ordered power : : can 
be expanded as 

(: X^ :) = J2 (t) (£'(-)'£'W'=-')e-*('=-20^_ (12) 

Hence, the moments can be obtained di- 

rectly by Fourier transforming the BHC data: 

^£(-)»^(+)m\^ _ y F("+™)(v3)e-*("-™)'^d^. (13) 
The moments of the original field are given by 

^£(-)n£(+)™^ = ^£(^)"£(+)™^e*^"~"'''"*. (14) 

Once we know how to measure the moments ^ in 
a single space-time point, the method can be extended 



to the general space-time correlations (Q. All field op- 
erators in the same space-time points can be grouped 
together and after a proper permutation any correlation 
function can be represented in the form 

a - ^ S n (r„i,)£W™nr.,t.) °o ^ , (15) 

where some of the rii or rrii can be zero. For example, 
the correlation function 

g^''^Hx,y,y)^{°o£^-\x)£^+^\y)°) (16) 

can be written in the form (|15|l as 

f;(i'^)(a;,y,y) = (S£(-)i(x)£(+)"(:r)£W"(z/)£W2(y)g). 

(17) 

Figure 13 illustrates the scheme for measuring the gen- 
eral correlation function H15|l . The phases ipi, . . . , ip^ 
of the local oscillator in each channel can be controlled 
independently by the corresponding phase shifters. For 
Tj — Vj^i — r and tj ^ tj^i we split the channel corre- 
sponding to the space point r into two parts as shown in 
Fig. ^ If more than two space points coincide, we split 
the signal in the corresponding channel by using the same 
tree-like combination of beam splitters and photodetec- 
tors as in MDd for a proper d, see Fig. ^ 

Measuring correlations of photodetectors of all the de- 
vices MD|^*'', one can extract the functions (|15|l with 
Ui + nii ^ 2^'^^, i = 1,...,N. Let us select ki pho- 
todetectors from the z-th device MD'-J^^ for i = 1, . . . , TV 
and measure the correlation function of all the K — 
ki + . . . + kN photodetectors chosen. One can easily see 
that such a correlation function F = p(fci - fciv) j-^^^g ^j^g 
form 

(18) 
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where P = J2i ^i^i and U is the number of those pho- 
todetectors of MD^'.-' that belong to the left subdevice 
Generalizing the single-mode approach, one can 
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compose the recorded data in the form 
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+ In, which is equal to 

N 

■ 1 



(20) 



Using multidimensional Fourier transform one 
can extract from the function p{ki...kN} _ 
^(fei...feiv)^^^^ _ (^^^ the correlations (|15|l or, equiva- 
lently, the original form of the space-time dependent 
field correlation function. 

Let us define the characteristic function C[u] via 



C[n] = 



exp 



where u — (ui, . 
derived by 



(21) 

.,un)- The moments (jJl are readily 



= (°£(-)"(a;,)fW'"(:E,)S>- (22) 



The general correlation functions (|15|l are given by the 
general partial derivatives of C[u] with respect to all the 
variables. Expanding the characteristic function as 
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it is seen that, in principle, the correlation functions pSf) 
completely describe the quantum statistical properties of 
the radiation field in the chosen space-time points. 

Finally, we note that our approach works only for 
quasimonochromatic fields that are completely recorded 
by the detectors. More generally, broadband fields could 
be studied by including spectral correlation measure- 
ments as was done for light pulses ■ This experiment 
also shows the feasibility of multichannel measurements 
as needed for our method. 

In conclusion, we have proposed a method for measur- 
ing general space-time dependent correlation functions of 
quantized radiation fields. Homodyne correlation mea- 
surements are performed and the data are combined and 



analyzed in a balanced form. The detected correlation 
functions are insensitive to imperfect detection. By using 
a strong local oscillator, even higher-order correlations 
can be determined. This opens new perspectives for the 
study of nonclassical correlations and, in particular, of 
entanglement of complex radiation fields. 
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